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Abstract 


This paper analyzes the asymptotic properties for small Bond number B of the equilib- 
rium capillary interface interior to a circular cylindrical tube vertically dipped in an 
infinite reservoir of liquid. (The Bond number B is a dimensionless parameter which is 
the ratio of gravitational to capillary forces.) The formal expansion in powers of B of 
the solution to the differential equation describing the equilibrium surf a^ (as can be ob- 
tained by standard perturbation methods) is proved to be truly asymptotic — to all orders 
and uniformly in the variable and parameter y, the contact angle. 

Sequences of general estimates, in closed form, from above and from below, are also 
given for the solution and related functions. The m^ term in these sequences are asymp- 
totically exact to order m. An idiosyncrasy of the problem, crucial in obtaining these 
estimates, is the absolute monotonicity of the structural function of the system in integral 
form. 



Figure 1. 
Introduction 


We consider the classical capillary problem of describing the equilibrium fluid inter- 
face 0 ^ interior to a circular cylindrical tube vertically dipped into an infinite reser- 
voir of liquid (Fig. 1) . Let u(r) be the height (above the level in the outer reservoir) 
of the surface as a function of the distance r to the axis of the tube. Then u(r) 
is a solution of the follo\;ing boundary value problem 


l(. - BU 

(1) 



u^(0+) = 0 

(2) 

u^(l-) « tan^j - vj. 

(3) 

this formulation, the quantities r, u, y, and B are dimensionless; 

Y is the con- 


2 

tact angle of s/ with the boundary cylinder (0 s Y * # and B ■ pgR /T > 0 is the Bond 

number with p the density difference across g the gravitational acceleration, R 

the radius of the cylinder, and T the surface tension. We refer for background to previ- 
ous papers and, in particular to^'^, for detailed proofs of most of the theorems. 

Boundary value problem (1,2,3) has a unique solution; this solution will be studied indi- 
rectly by the shooting method. Set v(r,h,B) to be the unique solution of the following 
initial value problem 


r 



2h + Bv 


( 4 ) 
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9r 


ii 


v(0,h,B) ■ 
v^(0»h,B) 


(5) 

(6) 


mean 


V 

Set aleo o(r,h,P) ■ — ■■ (r,h,B) and s * ain(j-Y)- (The parameter h le the 

curvature at the apex of the corresponding surface of revolution and o is the sine of the 
slope angle of v ) Then there exists a unique h ■ h(B»s) such that o(l,h»B) « s; this 
is also the unique h such that the function 2h/B v(-,hrB) is the solution of (1,2,3). 
Consequently, set 


u(r,B,a) - + v(r,h(B,a) ,B) . 


(7) 


For the above existence and uniqueness results, see, ^or exoinple, [6] or [2]. 

Instead of the initial value problem (4,5,6), we shall use an equivalent integral system 
for the pair (o,v) , namely, 


B <* 

o(r,h,B) ■ hr + - J pv(p,h,B)dp 
^ 0 

r 

v(r,h,B) * / g(o(p,h,B) )dp 
0 

where the structural function g is given by 


g(y) 




2n+l 


n-O 2 


for |y 1 < 1. 


( 8 ) 

(9) 

(10) 


Note that g is absolutely monotonic on the interval (0,1), i.c., has nonnegative deriva- 
tives on all orders on nSTlT^i This fact is essential in obtaining the estimates below. 

The approach will consist in first obtaining asymptotic results (resp. estimates) for the 
solution o(r,h,B), v(r,h,B) of the initial value problem (in integral form), then deriving 
from these results similar ones for h(B,s) and then u(r,B,s) . 


Absolute monotonicity and analyticity of o and v 

Since (4) is singular at r » 0, we need to show that the pair (o,v) is locally ana- 
lytic at r « 0 and depends analytically on the parameters h and B. This can be done in 
a number of ways (cf., for example, Wente^M ; alternatively^, it can be observed that the 
'one dimensional solution** (c,v) fi.e., the profile curve of the capillary surface bet%#een 
two vertical par * llel plates) is a majorant for (o,v) . Since (B,v) is solution of a 
regular initial v*alue problem, we conclude, by the method of majorants, that 

Proposition 1 ; For every tj, 8 >0, there exists p > 0 such that the functions 
o(r,h,B) and v(r,h,B) are analytic in the domain |r| < p, |h| < t\, |b| < ^ in 

Now the absolute monotonicity of g on (0,1) yields that o(r,h,B) and v(r,h,B) are 
absolutely monotonic in all three variables for positive values of r, h, and B. This can 
be checked by power series substitution into system (8,9). Thus the triple power series 
expansion of o and v at (0,0,0) has nonnegative coefficients. This last fact together 
with Pringsheim's Theorem below implies that this expansion is convergent on the maximal 
interval 0 < r < p » p(h,B) where o(*,h,B) and v(*,h,B) can be continued as a solution 
of integral system (8,9). Since o(*h,B) and v(*,h,B) are monotonic and bounded, we ob- 
tain convergence up to the boundary of the disc of convergence. 

Proposition 2 : The power series expansions of o(*,h,B) and v(',h,B) at r * 0 con- 

verge absolutely and uniformly in the closed convex disc |r| s p(h,B). Moreover 
0 < p(h,B) < • and, at r ■ p(h,£) , o(r,h,B) « 1. 

Coro] lary 3 ; The triple power series exoansions of o and v at (r,h,B) ■ (0,0,0) 
converge absolutely and uniformly in the closed domain 

D - ((r,h,B) ( C^: |r| * p(|h|jB|). (11) 

Prinqsheim^s Theorem ; Let 
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f(z) 


( 12 ) 


and R the radius of convergence of power series (12). Suppose c i 0, for all n ^ 0. 
Then z = R is a singular point of f. 

Asymptotic expansions 

Theorem 4 : For each m i 0, the series expansions in powers of B of the fxinctions 

h(B|S) , o(r,h(BrS) ,B ) , v(rrh(B,s) ,B) , and u(r,B,s) are asymptotic to order m uniform- 
ly in r and s over the entire range Oir^l, -l^ssl, as B— ^0. 

Proof (outline) : By symmetry, we may restrict ourselves to s i 0. Recall that h(B,s) 

is the unique solution of 

0(1, h,B) = s. (13) 


Now o(l,h,0)=h o(l,*,*) is analytic at each point (h,0) where h is in the com- 

plex disc Ih| <1. Since Oj^(l,h,0) = 1, the io^licit function theorem implies that 

h(B,s) is analytic at each point i0,s) where s is in the complex disc |sl < 1. A com- 
pactness argument shows that, for each s^ with 0 < Sq < 1, there exists Bq > 0 such 

that h(B,s) is analytic in 1 b! ^ B^, Is] ^ s^ and 

«D 

h(B,s) = E b” 5 (s) where 5_(s) = U4) 

n=0 " ^ 8 b" 


This yields the asymptotic statement for h(B,s) away frcnn s = 1, i.e., contact angle 

Y - 0. 

The neighborhood of s = 1 requires a special treatment and the use of the parametric 
system for the profile curve (parametrized by arclength) . The corresponding function 

h(B,a) where a = j - y is shown to be analytic in a neighborhood of (0,j) by the im- 
plicit function theorem- Thus h(B,s) == h(B,arcsin s) . We conclude, since the function 
arcsine is continuous at s = 1. 

The other three functions are handled in a similar way. QED. 

In particular, exists. This settles in the negative the problem of 

possible nonuniformity as s>l (cf. [1]). 

In Tables 1, 2, 3, and 4 below, we give the first few coefficient-functions of each ex- 
pansion as can be obtained by standard perturbation techniques. Previously Laplace^, 
Poisson®, Rayleigh®, and Concus® computed formally these asymptotic expansions to various 
orders. We should also mention that Siegel^® had proved recently the first term of those 
expansions to be truly asymptotic. 

First we need 

o(r,h,B) = Z (Br^)"x (hr) (15) 

n=0 ” 


v(r,h,B) = r I (Br^)^p (hr) 


n=0 


(16) 


(where the series converge in D) and recurrence formulas for the functions and p 

n n 


X,(x) 




(X/...X ") (i) 


1 ^ ^ (A, •••! ) 




(17) 

(18) 

(19) 


346 



wIiOM t:h0 aaoonci summation in (1!)) is Cox' x » p « iv* 

X n X * li 


Table 1, FiraU few X., 


Xjj(x) •■ X 


^ ,..v _ 1 X+aA-x‘ .. 

* 1 '’'' " J 7 r=?!y ’■ 




X^ (x) M 


tj) 


1 /vs i !> i 13 . U /r3 1 1 n 

" L Vl-x" 




'.fab).e a . First few u , 


Mq(x) • i (Wl-N-) 


L-(‘<A**X'' 


(-U (x) « -1 •(■ 3 log iisfcC 


„ i- 1 3- 1 ' S , 1 




$nbM 3 u f i^ut few 


§0 « Xq 


^1/’" :"^! V V " 
§2 » z * ^2 


H * ^i^a ^1^1** “ I ^1^1 




I 


Table 4 , Expansions for the BVP 


c(r,h(Bt?) ,B) = sr + B[r^Xj^(sr) - Xj^(s)] 

+ B^trXj^(s) X' (s) - r^X^(s) X' (sr) + r^X^Csr) ~ rX2(s)] + — 

v(r ,h(B,s) ,B) = rp^isr) + B[r^Pj^(sr) - r^X^(s) (sr) ] 

t B^[r^X^(s) X^(s)p^(sr) - r^ X^ (s) (sr) - ^ r^X^ (s) pj (sr) 

« r^Xj^(s)p^(sr) + r^^^isr)} + -•- 
u(r,B,s) = — + v(r,h(B,s),B) 

Estimates 

The coefficient- functions X^ and are odd, analytic in the unit disc, and abso- 

lutely monotonic in the interval (0,1) . This yields the following lower estimates by 
truncation of the series (15) and (16) , for each m ^ 0, 


m j 

o(r,h,B) i Z (Br^)^X (hr) 
n=0 " 

v(r,h,B) i r E (Br'^)*'n (hr) 
n=0 " 

for (r.h.B) € D+ = D fl [0,®)^. 

Using the fact that the right hand sxde of (20) is increasing in h and setting 
h = h(B,s) and r = 1 in (20) , we get the upper estimate 

h(B,3) s h (B,s) 
m 

where h^(B,s) is the unique nonnegative solution of equation 
m 

I B^X (h) = s- 

/% n 


( 20 ) 


( 21 ) 


( 22 ) 


(23) 


Estimates on the other side are obtained by induction on m from integral system (8,9), 
using Picard's method of successive approximations and a careful estimation of each iterate. 
The proof uses strongly the absolute monotonicity of the function g. For each m ^ 0, we 
get the upper estimates 


c(r,h,B) 


S ™E^(Br^)”A (hr) + (Br^)'"X^(o(r,h,B) ) 
n=0 " “ 


(24) 


v(r.h,B) s r E (Br^)"u (hr) +r (Br‘‘) M„(o(r,h,B) ) (25) 

n=0 " "* 

for (r,h,B) 6 d"'’. Using the fact that the right hand side of (24) is increasing in h and 
setting h = h(B,s) and r = 1 in (24) , we get, for mil, the lower estimate 

h(B,s) i h (B,s) (26) 

m 

where h (B,s) is 
m 

(i) 0, if B'"is/X^(s), 

m 

(ii) the unique nonnegative solution of equation 
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•* ■)« 


m-1 

1' B X (h) + B X (s) = s, 
n=0 " *" 

if b"* i s/X (s) . 

m 

Now, settinq h = h(B,s) in (20) and (21) and usinq (2t») , 

^ n 

o(r,h(B,s) ,B) i Z (Br") X (rh (B,sM 
n=0 " 


m 

v(r h(B,s),B) s r I (Br‘)‘V (rft (B,s)). 

n=0 "> 

Siiailarly, setting h = h(B,s) in (24) and (25) and using (22) , 

m~ 1 ^ ^ 

o(r,h(B,s) ,B) i I (Bi- )"X„(t-^ (B,s) ) + (Br‘)"'X^(sr) 
n=0 " ^ 

v(r.h(B,s) ,B) s. r'”r.\Br^)"M„(rh„ (B,s)) + r (Br^)"*u (sr) . 

n=0 ” "* 

Estimates (22 , 26 , 28 , 29 , 30 , 31 ) are valid for m:::!, B^O, O^s 
they are all asymptot ically exact to ordei m, as B 0. 


Some of these estimates had been obtained previously by Finn‘s 


s- 1 , and 0 2 ^ 
and Siegel^ ' . 


(27) 


(28) 


(29) 


(30) 


(31) 


r ^ 1; 
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